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Two-dimensional Problems in 
Polar Coordinates 



27 | General Equations in Polar Coordinates 

In discussing stresses in circular rings and disks, curved bars of narrow 
rectangular cros^ section with a circular axis, etc., it is advantageous to 
use polar coordinates. The position of a point in the middle plane of a 
plate is then defined by the distance from the origin 0 (Fig. 40) and by 
the angle 0 between -«• and a certain axis Ox fixed in the plane. 

/ Let us now consider the equilibrium of a small element 1234 cut out 
from the plate by the radial sections 04; 02, normal to the plate, and by 
two cylindrical surfaces 3, 1, normal to the plate. The normal stress 
component in the radial direction is denoted by the normal compo- 
nent in the circumferential direction by cr$, and the shearing-stress com- 
ponent by T r e, each symbol representing stress at the point r, 0, which is 
the midpoint P of the element. On account of the variation of stress 
the values at the midpoints of the sides 1, 2, 3, 4 are not quite the same 
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as the values a r , a e , r r e, and are denoted by (o>) x, etc., in Fig. 40. The 
radii of the sides 3, 1 are denoted by r 8 , r x . The radial force on the 
side 1 is (ov)i7i dd which may be written (ayr) i dd, and similarly the radial 
force on side 3 is — (o>r) 3 dd. The normal force on side 2 has a com- 
ponent along the radius through P of — (<r#) 2 (ri - r 8 ) sin (dd/ 2), which 
may be replaced by — (0-0)2 dr (dd/ 2 ). The corresponding component 
from side 4 is — (ere) a dr (dd/ 2). The shearing forces on sides 2. and 4 
give [(t t9 ) 2 — (T re )i\ dr. 

Summing up forces in the radial direction, including body force R per 
unit volume in the radial direction, we obtain the equation of equilibrium 

(<r r r)i dd — (oyr) 3 dd — (a e ) 2 dr~— (0-0)4 dr ~ 

+ [(tvs) 2 — (T r e)i ] dr + Rr dd dr = 0 

Dividing by dr dd this becomes 

(oVT") 1 (<r, -r) 3 1 r / \ 1 /- 'k i 1 (Jts) 2 (r r e) 4 , 77 ^ 

^ 2 H : — — + Rr = 0 

If the dimensions of the element are now taken smaller and smaller, to 
the limit zero, the first term of this equation is in the limit d(a T r) /dr. 
The second becomes oy, and the third dr r e/dd. The equation of equilib- 
rium in the tangential direction may be derived in the same manner. 
The two equations take the final form 



dov ! 1 dr r s , ov — <r 8 , D 
dr ^ r dd ^ r U 

1 dag dr r g . 2 t t o . a 

- -yv -) — ^ r r 

r dd dr r 



where S is the component of body force (per unit volume) in the tan- 
gential direction (^-increasing). 

These equations take the place of Eqs. (18) when we solve two-dimen- 
sional problems by means of polar coordinates. When the body force is 
zero they are satisfied by putting 

_ id# 1 d 2 <f> 

Cr r dr r 2 dd 2 

d 2 <j) 

= < 38) 
= i _ l = _ A A 

tB r 2 dd r dr dd dr \r dd J 

where 0 is the stress function as a function of r and 6. This of course 
may be verified by direct substitution. A derivation of (38) is included 
in what follows. 
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Instead of deriving (37) and observing that when R = S = 0 they are 
satisfied by (38), we can consider the stress distribution in question as 
first given in xy components <r x , <r y , r xy , as in. Chap. 3. We can then 
obtain from these the polar components o>, o>, T r g. From (13) we have 
(identifying a with B) 

cr r = cr x cos 2 B + a y sin 2 B + 2 r xy sin B cos 8 

o-g — cr x sin 2 B + <r y cos 2 B — 2r xy sin B cos B , (a) 

T r e = (cr y — a x ) sin 8 cos 8 + r^cos 2 8 — sin 2 8) 

We can similarly express a x , <x y , r xy in terms of o>, < r®, r rS by the relations 
(see Prob. 1, page 144) 

cr x = <x r cos 2 8 -f- o-g sin 2 8.— 2r r g sin 8 cos 8 

cr y = <j r sin 2 8 + eg cos 2 8 + 2r r e sin 6 cos 8 ( b ) 

T xy — (cr r — erg) sin 8 cos 8 -f- r r g( cos 2 8 — sin 2 8) 

To obtain (38) we consider next the relations between derivatives in the 

two coordinate systems. First we have 

r 2 = x 2 + y 2 8 = arctan - 

x 

which yield • - 



dr x dr y . n 

ox r x dy r 

dd V _ _ sin 8 dd _ x _ cos 8 

dx r 2 r dy r 2 r 

Thus for any function f(x,y), in polar coordinates f(r cos 8 , r sin 8 ), we 
have 

== 4 . ^ 1.^1 - df _ sin 8 df . 

dx dr dx "*■ dd dx C0S dr r dd ^ 

To get d 2 f/dx 2 we repeat the operation indicated in the last member of (c). 
Then 



cos 8 



sin 9 d s 
r dd . 



cos 8 



df sin 8 df N 



r dd j 



= cos 2 8 — — 5 — cos 8 sin 8 — 
dr 2 dr 



(~ a A 

\rddj 



sin 8 d f df\ sin 9 d ( . n df' 



With a little rearrangement, this takes the form 



d 2 f d 2 f 

= COS 2 8 -A + sm 



in 2 d(- ^ + i \ - 2 sin 9 cos 8 £ (- 

\r dr r 2 dd 2 / dr\rd8/ 
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Similarly, we find 






d 2 f 
dx dy 



sm 



y r Q r r 2 QQ2J 

e cos e (l a l + ± d -l - 

yr dr r 2 dd 2 dr 2 ) 



+ 2 sin 6 cos 6 — (- 
dr 



n*i) 

\rdd) 



- (cos 2 6 - sin 2 0) l (l §0 



He) 



(/) 



When we take for / the stress-function <t>(x,y) as in (29) — but with 
pg = o — the derivatives on the left-hand sides of (d), ( e ), and (/) become 
o-y, <?x, r xy , respectively. The expressions on the right-hand sides 'of (d), 
(e), and (/) can therefore be substituted for these stress components in the 
right-hand sides of (a). It is easily verified that the results reduce to (38). 

To convert the differential equation (a), page 35, to polar form, we 
first add ( d ) and ( e ) above to obtain 




dy 2 ) \dr 2 r dr r 2 dd 2 ) ^ 



(g) 



showing that the operator on the right is the polar equivalent of the 
laplacian operator on the left. Next, we find by addition of the first two 
of equations (b) 

<r* + CTy — ov + <70 Qi) 



For zero body force we have, as on page 30, 





+ <7y) ~ 0 



(0 



In view of (i), ( h ), and (g), this becomes 



> , 1 A , (<H , 1 dcf> , 1 d 2 <fi\ 

dr 2 r dr r 2 dd 2 ) \dr 2 ^ r dr ^ r 2 dd 2 ) U 



(39) 



From various solutions of this partial differential equation we obtain 
solutions of two-dimensional problems in polar coordinates for various 
boundary conditions. Several examples of such problems will be dis- 
cussed in this chapter. 



28 Stress Distribution Symmetrical about an Axis 



When the stress function 
bility (39) becomes 

AP + 1 d\ , 1 d± 

\dr 2 r dr) \dr 2 r dr 



depends on r only, the equation of compati- 



) 



d 4 <j> 2 d 3 4> 

dr 4 r dr 3 



i d 2 4> 
r 2 dr 2 



= 

r 3 dr 



0 (40) 
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This is an ordinary differential equation, which can be reduced to a linear 
differential equation with constant coefficients by introducing a new vari- 
) able t such that r = eK In this manner the general solution of Eq. (40) 
can be easily obtained. This solution has four constants of integration, 
which must be determined from the boundary conditions. By substi- 
tution it can be checked that 

<f> = A log r + Br 2 log r + Cr 2 -f- D (41) 

is the general solution. The solutions of a group of problems of sym? 
metrical stress distribution 1 with no body forces can be obtained from 
this. The corresponding stress components from Eqs. (38) are 

°> = 7 Ijr = 72 + B (! + 2 Io g r) + 2 C 
d 2 d> A 

(Te = dr* = ~ + 5 ( 3 + 2 lo S r ) + 2C (42) 

T r e = 0 

If there is no hole at the origin of coordinates, constants A and B vanish, 
since otherwise the stress components (42) become infinite when r = 0. 
Hence, for a plate without a hole at the origin and with no body forces, 
only one case of stress distribution symmetrical with respect to the axis 
may exist, namely, that when a r = ar g = constant and the plate is in a 
condition of uniform tension or uniform compression in all directions in 
its plane. 

If there is a hole at the origin, other solutions than uniform tension or 
compression can be derived from expressions (42). Taking B as zero , 2 
for instance, Eqs. (42) become 

- p + 2 C 

A (43) 

= - 7 + 2C 

This solution may be adapted to represent the stress distribution in a 
hollow cylinder submitted to uniform pressure on the inner and outer 
surfaces 3 (Fig. 41). Let a and b denote the inner and outer radii of the 
cylinder, and pi and p 0 the uniform internal and external pressur.es. 
Then the boundary conditions are 

(p’r)r—a Pi (p’ r ') r —b ~ Po ((I) 

1 The stress function independent of Q does not give all stress distributions inde- 
pendent of d. The function of the form Ad as in (g)‘on p. 126 illustrates this. 

2 Proof that B must be zero requires consideration of displacements. See p. 78. 

3 The solution of this problem is due to Lame, “Legons sur la theorie . . . de 
l’elasticit6,” Gauthier-Villars, Paris, 1852. 
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I 

■> 




Fig. 41 



Substituting in the first of Eqs. (43), we obtain the following equations 
to determine A and C: 

± + 2C=- Pi 
p+2 C=- V , 

from which 

A = a2h KVo ~ Vi) 
b 2 — a 2 

2 c = ~ Pob 2 

b 2 — a 2 




The radial displacement u is easily found since here ee = u/r, and for 
plane stress 



Ee 



6 ’ O’ Q VGf 



It is interesting to note that the sum o> + erg is constant through the 
thickness of the wall of the cylinder. ' Hence the stresses a r and a$ pro- 
duce a uniform extension or contraction in the direction of the axis of 
the cylinder, and cross sections perpendicular to this axis remain plane. 
Hence the deformation produced by the stresses (44) in an element of 
the cylinder cut out by two adjacent cross sections does not interfere 
with the deformation of the neighboring elements, and it is justifiable to 
consider the element in the condition of plane stress as we did in the. 
above discussion. 
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In the particular case when p 0 = 0 and the cylinder is submitted to 
internal pressure only, Eqs. (44) give 




(45) 



These equations show that o> is always a compressive stress and a 
tensile stress. The latter is greatest at the inner surface of the cylinder, 
where 






Pi(a* + b 2 ) 
b 2 — a 2 



(46) 



(o r «)max is always numerically greater than the internal pressure and 
approaches this quantity as b increases, so that it can never be reduced 
below pi, however much material is added on the outside. Various appli- 
cations of Eqs. (45) and (46) in machine design are usually discussed in 
elementary books on the strength of materials. 1 

The corresponding problem for a cylinder with an eccentric bore was 
solved by G. B. Jeffery. 2 If the radius of the bore is a and that of the 
external surface b, and if the distance between their centers is e, the 
maximum stress, when the cylinder is under an internal pressure pi, is the 
tangential stress at the internal surface at the thinnest part, if e < %a‘; 
and is of the magnitude 



Pi 



2 b 2 (b 2 + a 2 — 2 ae — e 2 ) 



(a 2 + b 2 ) (b 2 - a? - 2 ae - e 2 ) 
If e = 0, this coincides with Eq. (46). 



- 1 



29 | Pure Bending of Curved Bars 

Let us consider a curved bar with a constant narrow rectangular cross 
section 3 and a circular axis bent in the plane of curvature by couples M 
applied at the ends (Fig. 42). The bending moment in t his case is con- 
stant along the length of the bar and it is natural to expect that the stress 
distribution is the same in all radial cross sections, and that the solution 
of the problem can therefore be obtained by using expression (41). 

1 See, for instance, S. Timoshenko, "Strength of Materials," 3d ed., vol. 2, chap. 6, 
D - Va B- Nostrand Company, Inc., Princeton, N.J., 1956. 

2 Trans. Roy. Soc. (London), ser. A, vol. 221, p. 265, 1921. See also Brit. Assoc. 
Advan. Sci. Rept., 1921. A- complete solution by a different method is given in Art. 
66 of the present book. 

8 From the general discussion of the two-dimensional problem, Art. 16, it follows 
that the solution obtained below for the stress holds also for plane strain. 
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Fig. 42 



Denoting by a and b the inner and the outer radii of the boundary and 
taking the width of the rectangular cross section as unity, the boundary 
conditions are 



(1) 


CT r = 0 


(2) 


J a ere dr = 0 


(3) 


r r e = 0 



for r — a and r = b 
J a <rgr dr = —M 
at the boundary 



(a) 



Condition (1) means that the convex and concave boundaries of the 
bar are free from normal forces; condition (2) indicates that the normal 
stresses at the ends give rise to the couple M only, and condition (3) 
indicates that there are no tangential forces applied at the boundary. 
Using the first of Eqs. (42) with (1) of the boundary conditions (a) we 
obtain 



+ 5(1 + 2 log a) + 2(7 = 0 



A 
a 2 

^2 5(1 + 2 log b) -f- 2(7 = 0 



(b) 



Condition (2) in (a) is now necessarily satisfied. The use of a stress 
function guarantees equilibrium. A nonzero force-resultant on each end 
would violate equilibrium. To have the bending couple equal to M, the 
condition 



f b , [ b d 2 4> j 

/ ( Ter dr = / — — r dr = 
Ja Ja dr 2 

must be fulfilled. We have 



■M 



(d) 





64> b 
— — • r — 


[ b — dr = 


dcf> 

■ — • r 


Ja dr 2 


dr a 


Ja dr 


dr 
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and noting that on account of ( b ), 

d(f> b n 
— r =0 
dr a 

we find from ( d ), 

M = M 

or substituting expression (41) for <j>, 

A log ^ + B(b 2 log b — a 2 log a) + C(b 2 — a 2 ) = M 

CL 



(e) 



This equation, together with the two Eqs. (6), completely determines the 
constants A, B, C, and we find 



A = 



a 2 b 2 log - 
N . & a 



B = 



2M 

N 



(b 2 — a 2 ) 



C = 



M 



(/) 



jy [& 2 — a 2 + 2(6 2 log b — a 2 log a)] 

where for simplicity we have put 

N = (b 2 - a 2 ) 2 - 4a 2 6 2 ^log (g) 

Substituting the values (/) of the constants into the expressions (42) for 
the stress components, we find 

4:M l a 2 b 2 , & r . a\ 

r = ~ ~w \1^ log a + 6 log 5 + a log r) 



<re 
T r $ — 0 



- “ ^ (“ 7T lo S \ + &2 ] °g I + a 2 log ® + b 2 - a^j 



(47) 



This gives the stress distribution satisfying all the boundary conditions 1 
(a) for pure bending and represents the exact solution of the problem, 
provided the distribution of the normal forces at the ends is that given 
by the second of Eqs. (47). If the forces giving the bending couple M 
are distributed over the ends of the bar in some other manner, the stress 
distribution at the ends will be different from that of the solution (47). 
But, as Saint-Venant’s principle suggests, the deviations from solution 
(47) may be negligible away from the ends, say at distances greater than 
the depth of the bar. This is illustrated by Fig. 102. 

1 This solution is due to H. Golovin, Trans. Inst. Tech., St. Petersburg, 1881. 
The paper, published in Russian, remained unknown in other countries, and the same 
problem was solved later by M. C. RibiSre ( Compt . Rend., vol. 108, 1889, and vol. 132, 
1901) and by L. Prandtl. See A. Foppl, “Vorlesungen uber Technische Mechanik,” 
vol. 5, p. 72, 1907; also A. Timpe, Z. Math . Physik, vol. 52, p. 348, 1905. 
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m 




1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 
Values of ~r 



Fig. 43 



From the first of Eqs. (47) it can be shown that the stress <r r is always 
positive for the direction of bending shown in Fig. 42. The same can 
be concluded at once from the direction of stresses <t$ acting on the ele- 
ments n — n in Fig. 42. The corresponding tangential forces give 
resultants in the radial direction tending to separate longitudinal fibers 
and producing tensile stress in the radial direction. This stress increases 
toward the neutral surface and becomes a maximum near this surface. 
This maximum is always much smaller than (crg) max . For instance, for 
b/a = 1.3, (o>)max = 0.060 (cr*) max ; for 6/a = 2, (<r r ) max = 0.138M max ; for 
b/a = 3, (ov) max = 0.193(crj) max . In Fig. 43 the distribution of a e and 
<r r for b/a - 2 is given. From this figure we see that the point of max- 
imum stress <j t is somewhat displaced from the neutral axis in the direction 
of the center of curvature. 



30 | Strain Components in Polar Coordinates 

In considering the displacement in polar coordinates let us denote by u 
and y.the components of the displacement in the radial and tangential 
directions, respectively. If u is the radial displacement of the side ad 
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"of the element abed, (Fig. 44), the radial displacement of the side be is 
u + (du/dr) dr. The unit elongation of the element abed in the radial 
direction is then 






du 

dr 



(48). 



The strain in the tangential direction depends not only on the displace- 
ment v but also on the radial displacement u. Assuming, for instance, 
that the points a and d of the element abed (Fig. 44) have only the 
radial displacement u, the new length of the arc ad is (r + u ) dd and the 
tangential strain is therefore 

(r + u) dd — r dd _ u 
r dd r 



The difference in the tangential displacement of the sides ab and cd of 
the element abed is (dv/dd) dd, and the tangential strain due to the dis- 
placement v is accordingly dv/r dd. The total tangential strain is thus 1 



_ u . dv 
r r dd 



(49) 



Considering now the shearing strain, let a'b'c'd' be the position of the 
element abed after deformation (Fig. 44). The angle between the direc- 
tion ad and a'd' is due to the radial displacement u and is equal to du/r dd. 
In the same manner, the angle between a'b' and ab is equal to dv/dr. 
It should be noted that only part of this angle (shaded in the figure) 
contributes to the shearing strain and the other part, equal to v/r, 
represents the angular displacement due to rotation of the element abed 
as a rigid body about the axis through 0. Hence the total change in 
the angle dab, representing the shearing strain, is 



JrB = 



du j_ dV _ V 
r dd dr r 



( 50 ) 



1 The symbol eg was used with a different meaning in Art. 11. 



